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The polarization of the solar Mg ii h and k lines 
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^ ■ ABSTRACT 

O 

;_j , Although the h and k lines of Mg ii are expected to be of great interest for probing 

the upper solar chromosphere, relatively little is known about their polarization prop- 
erties which encode the information on the magnetic field. Here we report the first 
S::^ ! results of an investigation whose main goal is to understand the physical mechanisms 

that control the scattering polarization across these resonance lines and to achieve a 
|V] ! realistic radiative transfer modeling in the presence of arbitrary magnetic fields. We 

C/!) I show that the joint action of partial frequency redistribution (PRD) and quantum inter- 

'^ ! ference between the two excited y-levels produces a complex fractional linear polar- 

I ! ization (Q/I) pattern with large polarization amplitudes in the blue and red wings, and 

;-( ! a negative feature in the spectral region between the two lines. Another remarkable 

^ ! peculiarity of the Q/I profile is a conspicuous antisymmetric signal around the cen- 

ter of the h line, which cannot be obtained unless both PRD and 7-state interference 
effects are taken into account. In the core of the k line, PRD effects alone produce a 
triplet peak structure in the Q/I profile, whose modeling can be achieved also via the 
•/^ I two-level atom approximation. In addition to the Hanle effect in the core of the k line, 

t:j- I we emphasize also the diagnostic potential of the circular polarization produced by the 

Cf) ! Zeeman effect in the h and k lines, as well as in other Mg ii lines located in their wings. 
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1. Introduction 

It is through the Zeeman and Hanle effects that the magnetic fields of the solar atmosphere 
leave their fingerprints in the polarization of the emergent spectral line radiation (e.g., Casini & 
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Landi Degl' Innocent! 2007). The art of "measuring" solar magnetic fields thus relies on the de- 
velopment of suitable diagnostic tools based on such effects. The Zeeman effect produces circular 
and linear polarization signals whose amplitudes are proportional to "/? and 'J^, respectively, with 
H the ratio between the Zeeman splitting and the Doppler line- width. The Hanle effect is the mod- 
ification of the linear polarization produced by scattering processes in a spectral line, caused by 
the presence of a magnetic field. This line scattering polarization (whose amplitude is proportional 
to the fractional anisotropy of the incident radiation) is sensitive to magnetic strengths between 
approximately 0.25// and 55//, with Bh the critical Hanle field for which the Zeeman splitting 
of the line's level under consideration is equal to its natural width. Therefore, while the Zeeman 
effect in the far-UV (FUV) and extreme-UV (EUV) spectral lines that originate in the outer so- 
lar atmosphere is of limited practical interest for exploring its magnetism (because "/? «; 1), the 
Hanle effect is of great diagnostic potential because its magnetic sensitivity is independent of the 
wavelength and Doppler width of the spectral line under consideration. 

In this Letter we report the first results of a theoretical investigation aimed at clarifying the 
full diagnostic capabilities of the Mg ii h and k lines at 2802.7 A and 2795.5 A, respectively, 
with emphasis on the linear polarization pattern produced by scattering processes. These strong 
resonance lines are of interest for probing the hot-temperature plasma of the upper solar chromo- 
sphere (e.g., Uitenbroek 1997). The Mg ii h and k lines result from two transitions between a 
common lower level with angular momentum 7=1/2 and two upper levels with J = 1/2 and 



J = 311, respectively. As pointe d out in previous theoretical investigations (e.g. jAuer et al.lll980 : 



Belluzzi & Trujillo Buenoll201 1|) . quantum interference between such excited 7-levels is expected 
to produce observable effects in the wings of these lines, such as a cross-over of Q/I about the 
continuum polarization level in the spectral region between the two lines. This characteristic ob- 
servational signature of 7-state interference, pointed out by Stenflo (1980) in the Ca ii H and K 
lines, is expected to be more significant for the Mg ii h and k lines because their wavelength sep- 
aration and the ratio of the continuum to line opacity is m uch smaller than in the case of the Ca ii 



lines (e.g., see Fig. 5 of iBelluzzi & Trujillo Buendl201 ih . In fact, using the approximation of co 



herent scattering in the observer's frame, and simplifying the radiative transfer computations by 
assumin g that the line sour ce function and the integrated mean intensity are equal to the Planck 



function. lAuer et al.l (|1980|) estimated values of Qjl of about 4% in the positive polarization max- 
ima in the far wings, and of about - 1 .8% in the negative minimum between the two lines, for a line 
of sight (LOS) with jj. = cos 9 = OA (with 6 the heliocentric angle). 

Another key physical ingredient for modeling the Mg ii h and k lines is partial frequency 
redistribution (PRD). As shown by radiative transfer calculations carried out using a two-level 
model atom, PRD effects are expected to produce a conspicuous triplet peak Q/I signal in the core 
of the k line (e.g., Sampooma et al. 2010; see also figure 9 in Trujillo Bueno 2011, obtained in 
collaboration with Drs. M. Sampooma and J. Stepan). However, since the dispersion wings of 
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these lines are significantly overlapped and 7-state interference cannot be neglected, it is crucial to 
investigate the scattering polarization problem of the Mg ii h and k system accounting for the joint 
action of both PRD and 7-state interference effects. As we shall see below, our PRD with 7-state 
interference approach to this complex theoretical problem is not restricted to the line wings, so that 
our full radiative transfer calculations in model C of Fontenlaet al. (1993; hereafter FAL-C model) 
allow us to provide information on the amplitude and shape of the Q/I signal throughout the whole 
profile. A correct modeling of the scattering polarization in the core of spectral lines is particularly 
important because it is precisely in this spectral line region where the Hanle effect operates. We 
also provide estimates on the circular polarization induced by the longitudinal Zeeman effect. 



2. Formulation of the problem 



A quantum mechanical derivation of the redistribution matrix for polarized radiation, for a 
two-le vel a tom with unpolariz ed and infinitely sharp lower level, was carried out by lDomke & Hubeny 
(|l988h and lBommieii (|1997alJbb . In the atom rest frame (where the frequency and angular depen- 
dencies can be factorized), freque ncy redistributio n is described by the linear combination of Hum- 
mer's Rji and Rjjj functions (see lHummenll962l) . which describe purely coherent scattering and 
completely redistributed scattering, respectively. The corresponding branching ratios for the ^-th 
multipole component of the redistribution matrix are given by: 
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where Yr, F/ and F^ are the broadening widths of the upper level due to radiative decays, supere- 
lastic collisions (i.e., coUisional de-excitation), and elastic collisions, respectively, while D^^'' is 
the depolarizing rate due to elastic collisions (D^°^ = 0). Note that the branching ratios given in 
Eqs. (HI) and Q do not include the fa ctor (1 - e), with e = F//(Fr -i- F/) the photon destruction 
probability (cf., Eqs. (54) and (55) of Ib ommieii 1 1 997b|) . Their sum is equal to 1 for the ^ = 
multipole component (J3^'^^ = 1), while in the presence of depolarizing collisions it is generally 
smaller than 1 for K i^O. 

As mentioned in Section 1, quantum interference between different 7 levels is expected to 
produce significant observable effects on the scattering polarization signal across the Mg ii h and 
k lines. An atomic model acc ounting for the contribution of 7-state in terference (such as the two- 
term model atom presented in lLandi DeglTnnocenti & Landolfill2004i herefater LL04) is therefore 
needed for the investigation of this resonance doublet. Unfortunately, a rigorous quantum mechan- 
ical derivation of the redistribution matrix for such an atomic model, in the presence of elastic and 



inelastic collisions, has not been carried out yet. Nevertheless, the limit of coherent scattering in 
the atom rest frame is a reasonable approximation for modeling the scattering polarization across 
the Mg II h and k lines. This can be clearly seen from Fig. 1, where the branching ratio a corre- 
sponding to the k line (the one corresponding to the h line is practically identical) is shown as a 
function of height in the FAL-C model of the solar atmosphere. As shown in the figure, a is indeed 
equal to one at the heights where the line-core optical depth is unity for an observation at /z = 0.1, 
and it is very close to one also at the heights where the optical depth in the wavelength interval 
between the two lines is unity. 

A quantum mechanical derivation of the redistribution matrix for a two-term atom with un- 
polarized, infinitely sharp lower l evels, in the limit of purely cohe rent scattering in the atom's 



rest frame has been carried out by Landi DeglTnnocenti et al.l (|1997h within the framework of the 



metalevels theory. This redistribution matrix is suitable for investigating the scattering polarization 
across the Mg ii h and k lines, since the lower level of these lines is the ground level (it thus satisfies 
the hypothesis of an infinitely sharp lower level), and it has J{ = 1/2. Such redistribution matrix 
was obtained by neglecting collisi ons, and i s valid in the atom's rest frame. Following a derivation 
analogous to the one presented in iHummeii (|1962 r). we have calculated the corresponding expres- 



sion in the observer's frame, taki ng Doppler redistrib ution into account. An alternative derivation 



has been recently carried out by iSmitha et al.l (|201l|) starting from the Kramers-Heisenberg scat- 
tering formula. 

We consider a two-term model atom for Mg u, the lower term being composed of the ground 
level CS 1/2), the upper term of the upper levels of the h and k lines CP1/2 and ^^3/2, respectively). 
Although a rigorous derivation of the coUisional rates in a two-term atom has not been carried out 
yet, we included the eff"ect of (isotropic) inelastic and superelastic collisions in the redistribution 
matrix, under the following assumptions. We first assumed that the relaxation rate due to supere- 
lastic collisions is the same for the two upper 7- levels (see Section 7.13 of LL04 for the explicit 
expressions of such rate). We then assumed that such coUisional relaxation rate is also the one for 
the interference between the two upper 7- levels. Accordingly with the previous approximation, 
we calculated the inelastic collision transfer rates from the common lower level to the two upper 
levels. Since the lower level is unpolarized, these latter collisions only affect the populations of 
the upper 7-levels, and the corresponding rates are identical to the case of a multi-level atom (see 
Sect. 7.13 of LL04). CoUisional rates between pairs of 7-levels pertaining to the same term have 
been neglected, as well as elastic collisions. Under these hypotheses, following the convention 
according to which primed quantities refer to the incident photon, while unprimed quantities to 
the scattered photon, indicating with v and Q the photon's frequency (in the observer's frame) and 
propagation direction, we obtained 

^,/v.'.d';v,d) = 3§^J]J]J](-1)« 
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where P''^^(£)! , ^) . . is the ^-th multipole component of the scattering phase matrix (/, j = 0, 1 , 2, 3), 
Q is the scattering angle, Avo is the Doppler width (assumed to be the same for the two lines), and 
e' = TrIYj = AuclCut, with Auc the Einstein coefficient for spontaneous emission from the upper 
to the lower term, and Cut the superelastic collision relaxation rate. The function W is defined as 



W{a, v) = H{a, v) + iL(a, v) , 



(4) 



where H is the Voigt function, and L the associated dispersion profile. The reduced frequencies 
"J.J, and v'j^j^ are given by 

{yj„j, - v) , {vj^,j„ - y') 



"j.Jb - 



'JaJh 



(5) 



Avz) '''" Avo 

with Vj^ji^ the Bohr frequency between the levels J a and 7^. Although elastic collisions have been 
neglected (consistently with the assumption of purely coherent scattering in the atom rest frame), 
we took into account their broadening eff"ect in the evaluation of the broadening constant Y appear- 
ing in Eq. (|3]), and of the damping parameter a = Y/Avd- 



3. The scattering polarization pattern across the Mg 11 h and k lines 

In order to estimate, in the absence of magnetic fields, the amplitude and shape of the scatter- 
ing polarization pattern across the Mg 11 doublet, we have developed a non-LTE radiative transfer 
code based on the angle averaged expression (e.g. jRees & Saliballl982|) of the redistribution matrix 
described in the previous section, including the contribution of an unpolarized continuum. 

The numerical method of solution will be explained in detail in a forthcoming publication; it 
is based on a direct generalization to the PRD case of the Jacobian iterative scheme presented in 
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Trujillo Bueno & Manso Sainz (1999). The initialization of the iterative calculation was done using 
the self-consistent solution of the corresponding unpolarized problem, which we have obtained by 
applying Uitenbroek's (2001) radiative transfer code. This code has also been used to compute the 
inelastic and elastic coUisional rates, as well as the continuum total opacity (including the UV line 
haze contribution) and emissivity. 

In Figure 2 we show the results of the following PRD calculations of the 2/7 profile of the 
emergent radiation for a line of sight with ju = O.l: (I) the full two-term atom solution obtained by 
taking into account the impact of interference between the two excited 7-levels, (2) the two-term 
atom solution obtained by neglecting y-state interference effects, and (3) the solution for the Mg ii 
A;-line alone assuming a two-level atom with 7/ = 1/2 and 7„ = 3/2. 

The left panel of Fig. 2 shows the overall structure of the Q/I profile. Consider first the 
two-term atom solutions obtained by including (solid line) and neglecting (dashed line) 7-state 
interference. Clearly, the impact of 7-state interference is very important, since it leads to much 
larger polarization amplitudes in the wings (about 10% at around +15 A from the h and k line 
centers, respectively) and to a significant negative polarization feature between the h and k lines 
with a maximum amplitude of about -2%. The two-level atom solution for the k line (dotted line) 
produces wing polarization signals which are also significantly smaller than those corresponding 
to the full solution. 

The right panel of Fig. 2 shows the details of the Q/I profile in the core regions of the Mg ii 
h and k lines. The three solutions perfectly agree in the core of the k line, with a line-center 
amplitude of about 2%. However, it can be observed that the solution obtained by neglecting J- 
state interference (dashed line) does not show either the small asymmetry in the two negative Q/I 
peaks around the k-line core, or the clear antisymmetric Q/I feature that the combined action of 
PRD and 7-state interference effects produces around the center of the h line (whose polarizability 
is zero). 

Although the hypothesis of purely coherent scattering in the atom's rest frame is a good ap- 
proximation in the core of the lines and between them, frequency redistribution effects due to 
elastic collisions might not be completely negligible going towards the wings of the lines (see 
Fig. 1). In order to have a qualitative idea of the impact of this physical ingredient, we carried 
out a calculation including the contribution of the Rju redistribution function, under the limit of 
complete frequency redistribution (CRD) in the observer's frame. We derived the expression of 
Rjjj from the CRD theory presented in LL04, including the effect of inelastic and superelastic col- 
lisions as in the case of Rjj. We used the branching ratios a and (J3^^^ - a) defined in Eqs. (1) and 
(2) (strictly valid in the case of a two-level atom£|. We assumed D^^^ = (i.e. we neglected the 



'Note that the factor (1 - e) is included in our Rn and Rnj redistribution functions. 



depolarizing effect of elastic collisions). As expected, the contribution of Rju is negligible in the 
core of the lines, and very small between them. However, as shown in Fig. 3, it modifies in an 
appreciable way the amplitude of the Q/I profile in the wings of the lines. Note that the solution 
corresponding to the approximation of purely coherent scattering in the observer's frame is also 
shown in Fig. 3 (see the filled circles). 



4. Concluding comments 

As shown in this Letter, the physics that outside sunspots controls the scattering polarization 
of the Mg II h and k lines is the joint action of PRD and 7-state interference effects, which produces 
a complex Q/I pattern with sizable polarization maxima in the wings, and a negative polarization 
feature at wavelengths between the two line centers (see Fig. 2). These Q/I wing features can 
also be understood using the approximation of coherent scattering in the observer's frame (see the 
filled circles in Fig. 3o Another remarkable feature of the Q/I profile we have obtained is a clear 
antisymmetric Q/I signal around the center of the h line (see Fig. 2), which can only be found 
when taking into account the joint action of PRD and 7-state interference effects. In the core of the 
k line, PRD effects alone are responsible of the triplet peak structure of the Q/I profile (see Fig. 2). 

The reported linear polarization is sensitive to a magnetic field through the Hanle effect, which 
operates only in the core of the k line. Here magnetic fields weaker than 100 G are expected to 
produce significant changes in the line-center polarization amplitude. As shown in Fig. 2, it is a 
good new that the Q/I profile around the core of the Mg ii k-line can be well modeled through the 
two-level atom approximation, since this will facilitate the development of Hanle-effect diagnostic 
tools. In addition to the Hanle effect in the k-line, complementary information on the magnetic 
field of the solar chromosphere could be obtained by measuring the circular polarization produced 
by the Zeeman effect in all the Mg ii lines that are located between 2790 A and 2805 A (see caption 
of Fig. 1). The Stokes V signals of all these lines should be measurable. For example, a volume 
filling magnetic field of 100 G inclined by 45° with respect to the LOS produces V/I amplitudes of 
about 1% in all these Mg ii lines. 

Spectroscopy of the Mg ii h and k lines with n ovel space telescopes like the Interface Region 



Imaging Spectrograph (IRIS; see iTitle et al.ll2012|) will provide precious information on temper- 



atures, flows and waves. However, the magnetic field information is encoded in the spectral line 



^It is interesting to note that the observations performed by Henze & Stenflo (1987) could not confirm the negative 
polarization minimum between the h and k lines. Although the statistical significance of such pioneering observa- 
tions is very marginal, in a forthcoming investigation we will include the interaction with the continuum polarization 
processes in order to study its possible impact on the polarization amplitudes across these lines. 
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polarization, whose measurement with high spatial and/or tem poral resolution requ ires the devel 



opment of larger aperture telescopes such as SOLAR-C ( see iShimizu et alj 1201 11). ITie res ults 
shown in this Letter for the h and k lines of Mg n, and in iTrujillo Bueno et al.1 (|201ll |2012|) for 



the Lya lines of H i and He ii, strongly encourage the development of FUV and EUV spectro- 
polarimeters for the new generation of solar space telescopes. Additional results obtained in dif- 
ferent atmospheric models, along with the details of the methods used in this investigation, will be 
presented in a forthcoming paper. 

Financial support by the Spanish Ministry of Economy and Competitiveness through projects 
AYA2010-18029 (Solar Magnetism and Astrophysical Spectropolarimetry) and CONSOLIDER 
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Fig. 1. — The solid line shows the variation with height in the FAL-C model atmosphere of the 
branching ratio a of Eq. (1). The dotted line gives the atmospheric height where, for each wave- 
length, the ensuing optical depth is unity along a LOS with ju = 0.1, assuming that the only blends 
in the wings of the h and k lines are those produced by the following Mg ii lines: (a) two blended 
lines at 2797.930 A and 2797.998 A (resulting from the ^Pg/z - ^Dj/j and ^Pj^j - ^Dj/j transitions), 
and (b) one at 2790.777 A (caused by the ^Pi/2 - ^03/2 transition). Note that these absorption lines 
of Mg 11 originate a few hundred kilometers below the Mg 11 h and k lines, but they are still chro- 
mospheric lines. 
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Fig. 2. — The Qjl profile across the Mg ii h and k lines, calculated in the FAL-C model atmosphere 
for a LOS with ju = 0.1. Solid line: two-term atom PRD solution with 7-state interference. Dashed 
line: two-term atom PRD solution neglecting 7-state interference. Dotted line: two-level atom 
PRD solution for the fc-line transition. The reference direction for positive Q is the parallel to the 
nearest limb. The right panel shows in more detail the line core regions. 
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Fig. 3. — Comparison between the QII profiles obtained by considering only the Ru redistribution 
function (solid line, same as in Fig. 2), and a linear combination of Ru and Rm (dashed line). 
The filled circles give the Q/I values obtained under the approximation of coherent scattering in 
the observer's frame, after excluding the core regions where such an approximation to the full Ru 
solution does not hold. 



